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Abstract 



The complete density matrix for all polarization configurations in the process 
B* — ► ti, where B* is an off-shell Z or photon and t is the top quark, is calcu- 
lated numerically including oneloop QCD corrections, i.e. virtual and real gluon 
contributions in 0(a s ). The analysis is done in the framework of the helicity 
formalism. The results are particularly suited for top quark production at the 
Linear Collider, but may be useful in other circumstances as well. Relations to 
LEP and Tevatron physics are pointed out. 



1. Introduction 



Since its discovery in 1995 the top quark has been an object of increasing interest. 
The production process for top quarks has been analyzed in various theoretical 
studies both for pp and e + e~ collisions. Early references on the lowest order 
cross section are [0J for pp collsions and [[J for e + e~ annihilation. Higher order 
corrections to the cross section (total cross section, p T distribution etc.) have been 
calculated by several groups, [||, [|] for pp and ^ [J for e + e~ processes. These total 
cross sections do not involve information on the top quark polarization. Such spin 
effects only come in, if one studies distributions of top quark decay products. In 
some cases, spin effects have been studied, i.e. the distribution of the top quark 
spin vector, in ref. ]7j for pp and in ref. || for e + e~ collisions, but have not 
been extended to higher orders. An interesting step in this direction has been 
taken in ref. 0, where a Monte Carlo progam including final state spin terms 
has been written. That article is in fact concerned with higher order corrections 
to t + t~ production in e + e~ annihilation at lower energies and so does not take 
into account axial vector couplings. More general results, including axial vector 
couplings, can be found in ||, but not in the form of a Monte Carlo program. 
Ref. || is an alternative to the approach presented here, although we think that 
our approach is more enlightening, systematic and complete. 

Nothing is known about higher order spin corrections to the processes which in- 
duce top quark production in proton collisions (light quark annihilation qq — > ti 
and gluon-gluon fusion gg — > ti). In contrast to e + e~ annihilation, there is initial 
state gluon radiation and this makes those calculations very difficult. The simpli- 
fying feature of QCD corrections to heavy quark production in e + e~ annihilation 
is that the problem may be reduced to QCD corrections to the vectorboson- 
heavy quark vertex. Recently, a method for calculating QCD corrections to this 
vertex has been developed for the case of top quark decay t — > bW [[TIJ, i.e. the 
tbW vertex. In the present article this method will be applied and generalized 
to e + e~ — > B(= Z, 7)* — > ti. The method is based on the helicity formalism and 
allows to obtain the QCD corrections to the full spin density matrix of the pro- 
cess in a straightforward and economical way. Furthermore, the results obtained 
for top quark decay |TD] can be easily combined with the results of this paper by 



using a master formula Eq. (|]) to obtain QCD corrections for any distribution 
of top quark decay products in e + e~ annihilation one may be interested in. 

Within the Standard Model, all couplings of the top quark to other particles are 
completely fixed by its mass and by a few quantum numbers. For example, the 
coupling of the top quark to gluons is a pure vector coupling with strength g s , the 
coupling to a vector boson is given by ^7^ + 0^7^75 with v z = \ — az = \, 
w 7 = t^swCw, a 7 — etc. The couplings to the Z-boson and to the photon are 
particularly important for this article, because the process under consideration 
proceeds with an intermediate off-shell Z or photon where Z* and 7* arise from 
the annihilation of two massless fermions (e + and e~). As will be seen, this latter 
fact strongly reduces the number of independent helicity amplitudes and makes 
the results quite intuitive. 



2. Helicity Description 

Helicity amplitudes have been considered in many applications of phenomeno- 
logical importance in high energy physics, like jet production |Tl| , nonstandard 



effects in top quark processes | I2]| , and many others. The idea is, first to separate 



a given process into simpler subprocesses and then to explicitly evaluate all the 
possible spin amplitudes for the subprocesses in special Lorentz and Dirac frames. 
The results can afterwards be put together with the help of a master formula (to 
be given below). For example, consider the lowest order helicity amplitudes for 
top quark production in e + e~ annihilation through a vector boson B* with either 
vector or axialvector couplings ^ or to the top quark. They are given by 

M V (H, h, h) = VBUhipt^Vhipt)^ (1) 



M A (H, h, h) = a B u h {p t )l^v- h {pt)eH (2) 

where h = ±|, h = ±| and H = 0, ±1 label the spin states for the top quark and 
the B-boson. Note that for off-shell photons there is a longitudinal component 
with H = just as for a massive vector boson. In total, there are 24 amplitudes 
to be considered, which can in principle be relevant to the process. 



Real gluon emission cannot be treated on the amplitude level, but only on the 
level of the density matrix 

p u {H, h, h, H', h', hf) := M T (H, h, h)M}(H' \ hf, hf) (3) 

where I and J denote either V^(= vector) or A (=axialvector). The gluons are 
supposed to be unpolarized. After hermiticity one has 78 independent density 
matrix elements each for the VV, V A and AA total of 234 amplitudes. 

As will be shown below, only 60 of these are relevant for e + e~ annihilation. Note 
that the amplitudes are only determined up to an overall phase, and that this 
arbitrariness goes away when forming the density matrix elements. 

The full density matrix is in fact needed in the 'master formula', if one considers 
the combined production and decay process for the top quark. Assume the general 
case, that top quarks are produced in some process ab — > ti and then decay 
according to t — > W + b and i — > W~b, where the W's further decay to light 
(massless) fermions, W + — ► /1/2 and W~ — > 73/4. The cross section is then 
given by the 'master formula' 

0= IE M{h a ,h b ,h,h)M{h,h w+ )M{h,h w -)M{h w+ )M{h w -)\ 2 

EXT INT ,.v 

(4) 

where EXT = h a ,hf, denotes the spins of the external particles and INT = 
h, h, h w +, h w - the spins of the internal particles of the process. M(h a , h b) h, h) 
are the helicity amplitudes for the production process, M(h, h w -) for the decay 
of the antitop quark, and M{h w +) and M(hw~) are the amplitudes for the decay 
of the W + and W~ , respectively. The fact that the formula makes no explicit 
reference to the spins of the massless fermions fi, i=l,2,3,4, has the same reason 
as the non-appearance of the b-quark spins. Namely, the hf. are fixed by the 
V-A nature of the W decays, just as hf, is fixed by the V-A nature of the top 
decays (assuming massless b-quarks). Furthermore, if one is not interested in the 
decay of the antitop or of the W's, the corresponding amplitudes and helicities 
will not appear in the above formulas. In that case, the i and/or the W's will be 
one of the external (EXT) particles, whose spins have to be summed over after 
taking the square in Eq. ([|). Besides neglecting the b-quark mass, it is also 
a good approximation in Eq. (|J) to take the internal particles (top quarks and 
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Figure 1: Definition of the couplings 



W's) on-shell, because off-shell contributions are suppressed by powers of the 
width Y t and Yw 



3. Deconstruction of the e + e Production Cross Section 

In contrast, the intermediate vector boson Z and 7 must not be taken on-shell. 
The general spin amplitude for the production part of the process is 

M(h p , h e , h, h) = M^llhS + M,(h p ,K,h,h) (5) 

where h Vfi are the e spins and Q 2 is the square of the e + e~ energy. These 
amplitudes are now to be decomposed into the elementary amplitudes Eqs. ([!]) 
and (0). The decompostion will be on the algebraic level only, and no numerical 
analysis will be presented in this section, because we just want to show how the 
elementary amplitudes enter the e + e~ cross section. Afterwards the complete 
oneloop QCD corrections to all the elementary amplitudes resp. spin density 
matrix elements will be given (sections 5 and 6). 

The two amplitudes Mz and M 7 in Eq. (|J) have a decomposition of the form 

M B (h p ,h e ,h,h) = Y l M B (h p ,h e ,H) 1 2 M B (H,h,h) (6) 

H V ~ m B 

where the sum is over the Spins H = 0, ±1 of the off-shell B-boson. In this 
way the Bit vertex amplitudes Eqs. ([TJ) and (0) enter the cross section (Q). 



M B (h p , h e , H) are the spin amplitudes for e + e — > 5* and have a very simple 
structure due to the approximate masslessness of e ± . Namely, one has 



M B (h 



o: hei H) 



9 T2 for (h p ,h e ,H) = (+1-1,-1) 
§ for (h p ,h e ,H) = (-1+1+1) 
otherwise 



(7) 



where g± = v f ± af are the left- and right-handed part of the lepton couplings 



to the B, with i>f 



4 + 2s 



-- V 1 

2' e 



-2swCw, a 7 = 0. Furthermore, 



the spin quantization axis has been chosen to be the +z direction. According to 
Eq. (|7|) not all of the amplitudes M B (H, h, h) are picked up in Eq. (^) to give 
a nonzero contribution to the e + e~ amplitudes M B (h p ,h e ,h,h). Fig. 1 gives a 
graphic view about what notations are used for the couplings of the lowest order 
process. From this figure it becomes clear that the amplitudes Mb(H, h, h) are a 
linear combination ~ v B My + a B M A - Altogether one gets: 



M B (h p ,h e ,h,h) 



+^[v B M v 



l,h,h) + a B M A (-l,h,h)] h p = +\ 1 h 
-^[v B M v (+l, h, h) + a B M A (+l, h, h)] h p = - 



1 h 

2 ' 6 



otherwise 



2 

(8) 



One can now combine the contributions from 7 and Z (cf. Eq. ([5])): 



M(h p , h e , h, h) 



V-M v (-1, h, h) + a-M A (-l, h, h) h p = +±, h e = -~ 
v+M v (+l, h, h) + a+M A (+l, h, h) h p = -±, h e = +\ 
otherwise ^ 



where we have used the abbreviations v± = T^( g^t£p H — qt-) an d a ± 



1 9±az 



y/2Q 2 -m A 7 



According to Eqs. @ and (|9]), if the e beams are unpolarized, the 
cross section is a sum of two terms, a = o~_ + o~ + , where the first term is due to 
e + e _ helicities h p = +|, h e = —\ and the second term is due to e + e~ helicities 
h p = — |, h e = +|. One has 

a± = I ^[w±M y (±1, h, h) + o±M A (±l, h, h)]D(h, h)\ 2 (10) 



h,h 



where the contributions from the decay amplitudes have been summarized as 
D(h, ti). When one carries out the modulus squared in Eq. ([II]), it becomes ap- 
parent that in general all the density matrix elements of the form M/(+l, h, ti)Mj(+l, ti, h')* 
and M/(-l, h, h)Mj(-l, ti, ti)* for IJ = VV, VA + AV and AA are needed to 
calculate the cross section of the decay products. Note that instead of VA we 
are using the explicitly hermitean combination VA + AV. This is no restriction 
because to all orders in QCD one has MyM\ = MaM v . As will become explicit 
in section 5, the matrix elements of H = +1 and H = — 1 are related, so that 
only the first set M/(+l, h, h)Mj(+l, ti , ti)* has to be calculated. For each com- 
bination I J = VV, VA + AV and AA there are 10 independent matrix elements 
in this set. 

4. The Method 

The results to be presented were obtained in the helicity formalism. Furthermore, 
they concern the oneloop QCD corrections to the lowest order matrix elements. 
The lowest order expressions are usually simple, whereas the oneloop expressions, 
in particular for the case of hard gluons are quite lenghty for arbitrary spin orien- 
tations of B, t and i (the gluons are assumed to be unpolarized) []. Furthermore, 
the hard gluon contributions cannot be treated on the level of amplitudes, be- 
cause they have to be integrated over the gluon's energy and angles. One really 
has to go to the level of the (spin) density matrix, 

p(H, h, h, H', ti, ti) := M(H, h, h)M*(H', ti, ti) (11) 

to do the phase space integrations. There is, however, one circumstance which 
simplifies the task. This is related to the fact, that the oneloop QCD corrections 
to the total (spin-averaged) cross section for top quark production are known 
[gj, |14j| . This allows to get rid of the infrared and collinear singularities present 
in the matrix elements, by forming suitable singularity-free combinations of the 
spin-dependent and the spin-averaged expressions. The point is that the infrared 
and collinear singularities are 'universal', i.e. independent of the spin direction, 

^ee however ref. |]l3f for rather compact formulae 



so that they will drop out in suitable ratios and differences. To be more explicit, 
consider the ratio 

R(HM,^):=^ HM ^ K) (12) 

trp 

where the trace of the density matrix is given by 

trp = £ \M(H,h,h)\ 2 (13) 

H,h,h 

Now assume that each matrix element has the generic form p = T + aK where 
T stands for the tree level term, K for the higher order correction and a = CVfj 
is the QCD coupling constant and Cp = |. Assume similarly, that trp has the 
form trp = T a n + aK^. Then the ratio p/trp is given by 



R = = — + a KTan 9 TK ^ + 0{a 2 } = LO + a HO (14) 
trp Tf u 



It turns out that the difference KT & \\ — TK & \\ is free of infrared and collinear (and 
ultraviolet) singularities. If calculated, it gives an elegant means to determine 



higher order corrections to the density matrix [jig] . In the formulas and figures 
presented below, R will be written as R = LO + a HO where LO and HO will be 
plotted as a function of m t /Q. 

Finally note that the considerations in this section should in principle be done 
for the three matrices pvv, Pva(= Pav) an d Paa separately. In the case of V A it 
turns out that the trace of the density matrix vanishes (in leading order and in 
higher order). It is therefore convenient to consider a different linear combination, 
e.g. V A + AV + AA with a nonvanishing trp, and to apply the method Eq. (|Tjj) 
on that. Results for V A can afterwards be reconstructed by taking the difference 
of V A + AV + AA and AA (see Eq. (PT|) in the summary section). 



5. Lowest Order Amplitudes and Virtual Corrections 



I have calculated the eight needed amplitudes M/(+l, h, h), I — V, A as defined 
in Eqs. (0) and (Q) using the 'chiral representation' of 7 matrices, in which 



75 =diag(— 1, —1, 1, 1) etc. The results are presented in the rest frame of the 
virtual B-boson where the momenta are given by 

p B = (Q, 0,0,0) ft = |(l,0,0,/3) pf=|(1, 0,0,-/3) (15) 

with /3 2 = 1 — ^j-. An independent set of spinors for the top quarks is then given 
by 

u + i(pt) = (a+,0, a_,0) u_i(p t ) = (0, a_,0, a+) 

u+i(pf) = (-a+,0,a_,0) u_i(pr) = (0,-a_,0,a+) 

where a| = f (1 ± The quantization axis of the B-boson spin is chosen to 
be at an angle 9 w.r.t. the top quark momentum, i.e. the polarization vector for 
H — +1 is given by 

e+i = -4=(0, cos#, — i, — sin9) . (16) 
v2 

The most general parametrization of polarization vectors would be 

e_i = t=(0, cos cos # — « sin 0, sin (ft cos 9 + i cos0, — sin 0) 

v2 

e +i = _e -i 

e = — (0, sin# cos0, sin 6* sin 0, cos#) . (17) 

but for the problem at hand one may choose = without restriction. Note 
further that although the off-shell B-boson has a longitudinal polarization com- 
ponent, there is no contribution from longitudinal B's to the e + e~ cross section, 
as has been shown in section 3. Nevertheless, for completeness and curiosity, all 



the lowest order amplitudes are given here: 
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These amplitudes show a lot of symmetry. The most important for us arises from 
interchanging the role of particle and antiparticle. For = this amounts to 
9 <-» 6 + 71 and interchanges amplitudes of the form Mj(+1, h, h) <-> M/(— 1, /i, /t), 
I = V,A. In fact this relation is induced by CP invariance and it holds including 
higher order QCD corrections. It reduces the number of independent density 
matrix elements to be calculated by half from 60 to 30. 

Using Eq. ( |i8| ) one may calculate the trace of the corresponding density matrix 
Eq. (0): 



trp vv = 4 + 8— trpvA = trp AV = trp AA = A(3 



(19) 



These expressions should take the role of T a u in Eq. (|TJ). However, it appears 
that for the VA interference term the ratio Rva = Pva/^Pva becomes infinite. 



Instead of Rva the following combination will be considered in the following 

RvAAVAA '■= (PVA + PAV + PAa)/^PAA- 



The next step is to incorporate the corrections from virtual gluon exchange. This 
can be done on the amplitude level and is quite straightforward, because virtual 
gluons do not modify the kinematics of the lowest order process. The effect can be 
condensed to effectively change the vector and axialvector interactions according 
to 



7m -» 7m( x + a fi) + -^°v.vV V B a h 7m7s -> 7m7s(1 + af A ) 



(20) 



where a = as before and fx, f 2 and f A are functions of [3 and can be found 

in ||. It turns out that / 2 = fx — $a- Therefore, the amplitudes Eq. (|18|) can be 
easily extended to contain effects from virtual gluon exchange. We here give the 
result only for the relevant cases = and H — +1: 

M v (+1, -I -I) = sin 9(1 + af+) M A (+1, ~ ~) = 

M V (+1, -|, +|) = -^(1 + a/_) M A (+1, -I, +1) = -/jl±^(l + «/- 

M V (+1, +|, -|) = l^s4(l + a /_) M A (+1, +i -i) = + a /_ 



^sin#(l + a/ + ) M A (+l,+i,+i; 



M v (+1, +|, +|) = ^ sin 0(1 + a/+) M A (+1, +~ +1) = (21) 



where f + = f ± - £^f 2 and f_ = f A = f 1 - f 2 . 

It should be noted that the functions fx and f A merely 'renormalize' the form of 
the lowest order interactions, cf. Eq. (E^), and accordingly they do not contribute 
in the ratios Ru = pu/trpu, I J = VV,VAAVAA and AA. The interesting 
point is that the infrared and collinear singularities are solely contained in the 
functions fx and f A whereas f 2 is completely finite and given by 

1-02 1 _ R 

f- = ^ la r4 (22) 

The normalized density matrix elements i?yy(+l, h, h, +1, h', h!) are given in the 
table where A = 1-2m t/Q s = s [ n Q anc i c = cos Q anc l the *'s follow from the 

mt/Q ' 1 1 
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symmetry of the matrix. It is thereby explicit that all corrections are ~ f 2 as 
anticipated. 

The density matrix elements of Raa and Rvaavaa do not get any corrections 
at all from virtual gluon exchange because according to Eq. (^l|) there are two 
factors of 1 + ct/_ both in the numerator and denominator of Paa/^Paa and 
similarly for Rvaavaa- 



6. Real Gluon Emission and Numerical Results for the Corrections 
to the Normalized Density Matrix 

Next we come to the two Feynman diagrams with real gluon emission B* — > tig. 
The higher order corrections are calculated in such a way that all gluon d.o.f. 
are summed and integrated over. This means, for example, the gluon is assumed 
to be unpolarized. Furthermore, a rather complicated phase space integrattion 
has to be performed. In contrast to the tree level process B* — > ti where the 
phase space is trivial, one has here two highly nontrivial integrations which I 
have choosen to perform numerically. This is then straightforward because, as 
pointed out before, the integrand corresponding to Eq. is completely finite. 



The trivial phase space for the lowest order kinematics (|T5|) is given by PS(B* 



ti) — III ^7^{Pt + Pi — Pb) = f P- With an additional gluon the 4-momenta of 
the particles become more complicated. First, the energies of t and i are not fixed 
as in lowest order - in fact they are to become the integration variables - and 



secondly, there is now an angle uj between the top quark and antitop direction. 
Pb = (Q, 0, 0, 0) p t = ^zi(l,0,0,/3i) p? = ^:r 2 (l,/3 2 sinu;,0,/3 2 cosa;) 

(23) 

Here xi and x 2 are the (normalized) energies of t and i and (3f : = 1 — Jzqj, z = 1, 2. 
w is not an independent variable but can be related to x\ and x 2 through the 
gluon's on-shell condition p 2 g = and 4-momentum conservation p s = PB—Pt—pi- 
One obtains: 

TYl 2 I 

= 1 + - xi - x 2 + -x lX2 (l - P1P2 coscu) (24) 
The 2-dimensional phase space integral is given by 

PS(5* -> %) = — ^- / dy dz (25) 



4 Jo 



where 



2/ := — = 1 - Xi z := — = 1 - x 2 (26) 



and 



^ = ^Ml±A)-2^] (27) 



The polarization vectors Eqs. (|16D and ( |TTD of the B-boson are not modified 
in higher orders, because ps has not changed. In contrast, the spinors for the 
fermions become more complicated than in lowest order: 

u + i(pt) = (oi + ,0,ai_,0) u_i(p t ) = (0,ai_,0,ai+) 

«+|(Pf) = («2+S2, -fl2+c 2 , -a 2 _s 2 ,a 2 _c 2 ) u_|(pf) = (a 2 -C 2 , a 2 _s 2 , -a 2+ c 2 , -a 2 

where s 2 = sin ~, c 2 = cos ^ and a 2 ± = ± i = 1, 2. 

Using these parametrizations, the real gluon amplitudes Mi(B* — > tig)(H, h, h), 
I = V,A and h,h — ±|, must be calculated. As in lowest order, it is sufficient 
to know the case H — +1. Before performing the phase space integrations, the 
density matrix must be formed. More precisely, we have calculated the normalized 



density matrices Ru(+1, h, h, +1, ti, ti), IJ = VV, V AAV AA and AA, Eq. (TTJ), 
corresponding to the finite combination Eq. (|T4|). Numerical results for all the 
30 independent matrix elements are shown in Figures 2-61 as a function of the 
e + e~ energy Q and for a top quark mass of m t = 170 GeV. The results are shown 
in the following form: each density matrix element has a decomposition of the 
form 

R u {+1, h, h, +1, ti, h!) = (7^(1 + c?) + C 2 s 2 t + C 3 c t + C 4 s t + C 5 c t s t 

(28) 

where c% = cos 9, St = sin 8 and of course the coefficients Cj depend on the spin 
quantum numbers h, h, ti and ti. Each coefficient Cj has a lowest order con- 
tribution, which however vanishes in many cases, and a higher order correction, 
i.e. 

Cj = LO + a HO (29) 

where a = Cf^ as before. In Figures 2-61 there is always one figure displaying 
the lowest order term (LO) for each j = 1, 2, 3, 4, 5 followed by the corresponding 
figure for the higher order term (HO); and these pairs of figures are successively 
presented for each combination (+1, h, h, +1, ti, ti). Note that only the following 

combinations of (h, h, ti, ti) are shown: ( ), ( h), ( 1 — ), ( 

++), (- + -+), (- + +-), (+ - +-), (+ + -+), (+ + +-), (+ + ++). The rest 
follows from the symmetry of the density matrix, i.e. Ru(+1, h, h, +1, ti, ti) = 
Ru(+1, ti, ti, +1, h, ti). Some of the curves seem to increase rather strongly with 
increasing Q 2 . However, it can be shown that all HO corrections converge to a 
constant at Q — > oo. Corrections are typically at the percent level. This is not 
astonishing, because by forming the normalized density matrix, the bulk of the 
higher order correction due to the known overall K-f actors [TJJ] drops out. The 
odd figure numbers 3-61 really contain the main results of our paper, from which 
higher order QCD corrections to any angular distribution of top quark spins or 
top quark decay products can in principle be determined (for the latter one needs 
in addition the QCD corrections to the top quark decay amplitudes calculated in 

ED- 



For convenience, we present in figures 62-91 the 6 dependence of the matrix 



elements Ru(+1, h, h, +1, h', h'). More precisely, the matrix elements are written 
as Rjj(+1, h, h, + 1, h', h') = LO + a HO and LO and HO are displayed as a 
function of 9 for a fixed total energy Q = 400 GeV. Again, the corrections are 
relatively small because the normalized density matrix is considered, in which 
the correction to the trace of the density matrix drops out. 

7. Summary 

In this article a complete decomposition into spin contributions of the production 
of top quarks in e + e~ annihilation has been given including higher order QCD 
corrections. The corrections are relatively small (on the percent level) because 
the normalized density matrix was considered, in which the correction to the trace 
of the density matrix drops out. This means in turn that twoloop calculations 
are at most necessary for the trace of the density matrix (corrections of up to 
10% |14|), but not for any of the normalized density matrix elements, because 
these 0(a 2 s ) corrections are most probably below the permille level. 

Our results can be applied to calculate distributions for any decay product of 
top quarks at a future linear collider. I have not attempted to work out these 
applications, but have merely given in section 3 the necessary reduction formulae. 
To prove that the applications are rather straightforward, I would like to finish 
this work by making the QCD corrections to the e + e~ cross section Eq. (|i~0|) 
explicit. Eq. ( JlOj ) can be rewritten as 

a±= ]T D(h,h)D*{ti,h') 

h,h,h',h' 

2 

c — — TfL 

{vlR vv (±l, h, h, ±1, ti, h'){A + 8— 1)(1 + aK vv (f3)) + 4/? 2 (l + oKaaW) 
[v±a±R VAAVAA (±l } h, h, ±1, ti, h!) + a±(a± - v±)R AA (±l, h, h, ±1, ti, htffy) 

where KyviP) and K AA {j3) are the known corrections to the trace of the density 
matrix [[14]], cf. Eq. (|l^). Note that Ky A = K A y = 0, so that there is an overall 
factor of l + aK AA in the second line of Eq. ([TI]). Using Eq. fl3"T| ) one can directly 
insert the corrections (figures 2-61) to obtain the cross sections including higher 
order QCD. 
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